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Abstract 

We develop a novel method to identify active galactic nuclei (AGNs) and study the nature of low-luminosity 
AGNs in the Sloan Digital Sky Survey. This is the second part of a series of papers and we study the correlations 
between the AGN activities and host galaxy properties. Based on a sample of AGNs identified with the new method 
developed in Paper-I, we find that AGNs typically show extinction of ry = 1.2 and they exhibit a wide range of 
ionization levels. The latter finding motivates us to use [oil] + [oill] as an indicator of AGN power. We find that 
AGNs are preferentially located in massive, red, early-type galaxies. By carefully taking into account a selection bias 
of the Oxygen-excess method, we show that strong AGNs are located in actively star forming galaxies and rapidly 
growing super-massive black holes are located in rapidly growing galaxies, which clearly shows the co-evolution of 
super-massive black holes and the host galaxies. This is a surprising phenomenon given that the growths of black 
holes and host galaxies occur at very different physical scales. Interestingly, the AGN power does not strongly 
correlate with the host galaxy mass. It seems that mass works like a 'switch' to activate AGNs. The absence of 
AGNs in low-mass galaxies might be due the absence of super-massive black holes in those galaxies, but a dedicated 
observation of nuclear region of nearby low-mass galaxies would be necessarily to obtain deeper insights into it. 
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1. Introduction 

Massive black holes were first speculated as powerhouses 
for luminous quasars (e.g., Lynden-Bell 1969). A large amount 
of work conducted afterwards vastly improved our understand- 
ing of active galactic nuclei (AGNs) and it is now widely rec- 
ognized that super-massive black holes are far from rare, exotic 
objects and they are basic constituents of massive galaxies. A 
large spectroscopic survey of nearby galaxies has been a suc- 
cessful strategy to unveil statistical nature of AGNs. Ho et al. 
(1995) reported on a systematic survey of 486 nearby galaxies 
with the Hale 5m telescope and Ho et al. (1997) studied weak 
AGNs from the survey, most of which turned out to be low- 
ionization objects often called LINERs (Heckman 1980). They 
found that AGNs are preferentially located in massive early- 
type galaxies with large bulges. Further observations of nearby 
galaxies revealed that super-massive black hole mass correlates 
with host galaxy properties, particularly with bulge luminos- 
ity and mass (Kormendy & Richstone 1995; Magorrian et al. 
1998; Ferrarese & Merritt 2000; Gebhardt et al. 2000). This 
suggests that the super-massive black hole growth and galaxy 
growth are linked, and it has motivated a lot of effort to study 
the black hole evolution in the context of galaxy evolution both 
theoretically and observationally. 

Recent galaxy formation models, which incorporates energy 
feedback by AGNs in one way or another, seem to reproduce 
observed properties of galaxies reasonably well. AGNs seem to 
be a favorable energy feedback mechanism due to thek contin- 



uous energy output in massive galaxies with QSO-like strong 
feedback in starburst phase (e.g., Granato et al. 2004; Springel 
et al. 2005; Bower et al. 2006; Croton et al. 2006), which helps 
reproduce the 'break' in the massive end of the stellar mass 
function. However, these simulations cannot directly follow 
the AGN physics, which happens on less than a parsec scale, 
and rely on a simplified recipe. Also, the exact form of en- 
ergy feedback is still highly uncertain. Despite the simplicity 
and assumptions employed, however, these models succeeded 
to reproduce overall properties of observed galaxies relatively 
well. 

Observational evidence of AGN feedback, on the other hand, 
is still limited and we do not yet know how exactly AGN af- 
fects the galaxy evolution. High-quality large spectroscopic 
surveys of the local universe such as the 2dF survey (Colless 
et al. 2003) and the Sloan Digital Sky Survey (SDSS; York et 
al. 2000) rapidly improved our understanding of AGNs. It is 
now widely accepted that AGNs are preferentially located in 
massive galaxies (Kauffmann et al. 2003). Strong AGNs tend 
to reside in blue galaxies and high ionization AGNs (Seyferts) 
also slightly prefer blue galaxies over red galaxies (Kewley et 
al. 2006). It seems that very strong AGN activities may be af- 
fecting star formation activities of the host galaxies. Ho (2005) 
studied nearby QSOs and found possibly suppressed SFRs of 
the host galaxies. Kim et al. (2006) extended the analysis with 
a large number of type-I AGNs drawn from SDSS and sug- 
gested that the observed [oil] emission can be explained en- 
tirely by the photo-ionization due to AGN with little emission 
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from HII regions. Recently, Greene et al. (2011) examined 15 
luminous QSOs and reached essentially the same conclusion. 
If we turn our attention to low-luminosity AGNs, which com- 
prise the majority of the overall AGN population, some authors 
claimed that the host galaxies tend to show colors in between 
the red sequence and blue cloud (so-called green valley), which 
led them to speculate on on-going effects of AGN feedback 
(Schawinski et al. 2007; Schawinski et al. 2010). Heckman et 
al. (2004) studied the mass of the AGN host galaxies and sug- 
gested that most of the active super-massive black holes are 
located in intermediate mass galaxies, therefore representing 
down-sizing in AGN activities. 

A significant fraction of the previous studies from the sur- 
veys of the local universe is based on optical emission line 
diagnostics first proposed by Baldwin et al. (1981; here- 
after BPT) to identify AGNs thanks largely to a large number 
of high quality spectra from SDSS. However, a major draw- 
back of the popular BPT diagnostics is that it requires four 
emission lines which are not always easy to measure at high 
signal-to-noise ratios. Among them, H/3 can be fairly weak 
in AGNs and that imposes a practical limit on its sensitivity 
to low-luminosity AGNs (Paper-I). Several authors proposed 
variances of the BPT diagnostics to overcome the issue (Rola 
et al. 1997; Lamareille et al. 2004; Yan et al. 2006; Yan et al. 
201 1 ; Juneau et al. 201 1), but they may not be fully satisfactory 
alternatives. In particular, some of them make a priori assump- 
tion of host galaxy properties and therefore may not be suited 
to study AGN - host galaxy correlations. 

To improve the situation, we have developed a novel method 
to identify AGNs in Paper-I based on a very simple idea of 
comparing expected with observed emission line luminosities. 
We do not make any a priori assumptions of host galaxy prop- 
erties. While the method suffers from contamination and in- 
completeness as all the other AGN identification methods do, 
it is a statistically efficient method to identify low-luminosity 
AGNs. With this new method, we aim to study the nature of 
low-luminosity AGNs and their host galaxy properties in this 
paper 

The structure of this paper is as follows. In Section 2, we 
give a brief overview of the new AGN identification method 
developed in Paper-I. We then develop a basis to characterize 
AGN activities in Section 3, followed by a detailed study of 
correlations between AGN activities and host galaxy proper- 
ties in Section 4. We summarize our results and discuss their 
implications in Section 5. Unless otherwise stated, we adopt 
= 0.3, riA = 0.7, and Hq = 70 km s^^ Mpc^^. All the 
magnitudes are given in the AB system. We use the follow- 
ing abbreviations ; AGN for active galactic nucleus, BPT for 
the Baldwin et al. (1981) diagnostics, SF for star formation 
and SFR for star formation rate. Emission lines used in this 
work include [oil] AA3726,3729, UP A4861, [oiii] A5007, 
[oi] A6300, Ha A6563, [Nil] A6583, and [sii] A6716,6730. 

2. Overview of the new method and the data 

2.1. Method 

We briefly give an overview of the new method to identify 
AGNs. For the full development, readers are refereed to Paper- 
I. 
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The most widely adopted optical emission line diagnostics 
to identify AGNs is the one proposed by Baldwin et al. (1981) 
and later extended by Veilleux & Osterbrock (1987). It is based 
on the idea that an ionizing spectrum of AGN is often harder 
than that of young stars, and emission lines show characteristic 
intensity ratios. However, in principle, a single emission line 
contains information about AGN. An emission line luminosity 
of a galaxy we measure from its spectrum is 

Lmeasured = LsF + LaGNi (1) 

where L5 p and Lagn are luminosities due to star formation 
and AGN, respectively. It has been a long standing problem 
to separate the two components, which has hindered detailed 
studies of pure AGN emission. But, we have developed a 
novel method to estimate Lsf in Paper-I. AGNs often emit 
strong emission lines, but their continuum is generally much 
weaker. Optical continuum emission of a galaxy spectrum is 
dominated by stellar light, not by AGN, in most AGNs (Schmitt 
et al. 1999). The stellar continuum contains information about 
star formation rates (SFRs) and dust extinction of the galaxy. 
We extract these two pieces of information by fitting model 
templates generated with an updated version of the Bruzual 

6 Chariot (2003) code to the observed spectrum. From SFRs 
and dust extinction, we can work out Lsf- We then compare 
Lraeasured with LsF- If wc obscrvc a significant Lagn, a 
galaxy likely hosts an AGN. As we will discuss in Section 3, 
AGNs exhibit a wide range of ionization state and we use a 
sum of [on] and [oill] luminosities to identify AGNs. We 
have shown that we can estimate Lsf for [Oll]H-[oill] with a 
factor of ^ 1.7 accuracy in Paper-I. Since we identify AGNs 
using excess Oxygen emission, we dub the method "Oxygen- 
excess method". Note that we will extensively use stellar mass, 
SFR and dust extinction from the spectral fits in this paper. 

Compared to the commonly used Baldwin et al. (1981) di- 
agnostics, this Oxygen-excess method gives the same SF/AGN 
classifications for 85% of objects that exhibit strong enough 
emission lines to apply the BPT method. The BPT diagnostics 
uses 4 emission lines and it takes ratios of them, requiring a 
high signal-to-noise ratio for each line. On the other hand, the 
new method uses a sum of the two lines and its requirement for 
signal-to-noise is less demanding. As a result, the method is 
applicable to 78% of the objects in our sample, while BPT is 
applicable only to 43%. Our method is thus fairly sensitive to 
low-luminosity AGNs. We have quantified average properties 
of the Oxygen-excess objects by stacking the spectra and we 
have shown that many of the Oxygen-excess objects, if not all, 
likely host AGNs. Serendipitous X-ray observations (Evans et 
al. 2010) as well as radio data from the FIRST survey (Becker 
et al. 1995; White et al. 1997) also show that SF/AGN classi- 
fications are made well. Another unique feature of the method 
is that it allows us to subtract emission line luminosities due 
to star formation and extract AGN luminosities, which are cru- 
cial to characterize AGNs. We will make a full use of these 
unique features to study the nature of low-luminosity AGNs in 
this paper 

2.2. Data 

We use data from the Sloan Digital Sky Survey Data Release 

7 (Abazajian et al. 2009). Using a dedicated 2.5m telescope 
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(Gunn et al. 2006), the SDSS surveyed a quarter of the sky in 
both imaging and spectroscopy (Fukugita et al. 1996; Gunn et 
al. 1998; Eisenstein et al. 2001; Strauss et al. 2002; Richards et 
al. 2002; Gunn et al. 2006; Doi et al. 2010). The galaxy sample 
for our study is drawn from the main galaxy sample (Strauss et 
al. 2002), which is a flux-limited sample down to r = 17.77 se- 
lected from the imaging survey. We apply the following selec- 
tion criteria; SPECClass=2 (i.e., objects are galaxies) located 
at 0.02 < z < 0.10 with high confidence flags (zCONF>0.8 
and zWarning=0). We intentionally remove QSO-Uke ob- 
jects (specClass=3) because our method is not applicable to 
those objects (continuum needs to be dominated by stars, not 
by AGN), leaving 283,031 objects in total. We apply a first- 
order correction for the fiber loss to the stellar mass and SFR 
by assuming that light in the fiber is representative of the en- 
tire galaxy. The corrected quantities are shown with a subscript 
apercorr. 

2.3. Caveats on photo-ionization due to non-AGN sources 

Before we present our results, it is important to remind of the 
contamination of objects that are photo-ionized by non-AGN 
sources. As discussed in Paper-I, there are a number of pos- 
sible non-AGN sources that could produce LINER-like emis- 
sion line ratios. They include shocks (Cox 1972; Heckman 
1980; Dopita & Sutherland 1995; Dopita & Sutherland 1996), 
Wolf-Rayet stars (Terlevich & Melnick 1985; Kewley et al. 
2001), post-starburst (Taniguchi et al. 2000) and post-AGB 
stars (Binette et al. 1994; Stasiriska et al. 2008; Sarzi et al. 
2010; Cid Fernandes et al. 2011). Relative contributions of 
these ionizing sources to the observed LINER populations are 
uncertain, but they may contaminate our AGN sample. 

As shown later, most of the Oxygen-excess objects are red, 
massive galaxies. In these galaxies, shocks, Wolf-Rayet stars, 
and post-starburst are unlikely the primary origin of the ob- 
served emission. Post-AGB stars remain a concern. Several 
attempts have been made to quantify their role (Binette et al. 
1994; Stasiriska et al. 2008; Sarzi et al. 2010; Cid Fernandes 
et al. 2011; Yan & Blanton 2011). In particular, Sarzi et al. 
(2010) and Yan & Blanton (201 1) studied nearby galaxies and 
observed spatially extended Ha or H/3 emission, which is un- 
likely to be fully explained by a central source. There is no 
doubt that these post-AGB stars contaminate the true AGN 
population and the observed emission line luminosities of the 
Oxygen-excess objects are at least partly contaminated by post- 
AGB photo-ionization. 

We argue in Section 3.3 that emission due to post-AGB stars 
does not seem to dominate the overall emission, but the ex- 
action fraction of objects that are photo-ionized by post-AGB 
stars among the Oxygen-excess objects and the exact contri- 
bution of post-AGB stars to the observed emission remain un- 
clear As discussed in Paper-I, the SDSS data are not a right 
data set to precisely pin down the role of post-AGB stars and 
we do not try to pursue the issue further in this paper. But, we 
should bear in mind that this unknown role of post-AGB stars 
remains one of the major uncertainties throughout the paper. 
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Fig. 1. Distribution of tlie Balmer decrement. The star formation com- 
ponent is subtracted from the Ha and H/3 luminosities and extinction 
due to the interstellar medium is corrected for. The dashed lines show 
Ty assuming the dust-free Balmer decrement of 3.1. Note that one 
can convert ry into the commonly used extinction in the V-band using 
Ay = 1.09t\/. We use objects with Ha and H/3 detected at more than 
lOcr only here to quantify the line ratio. 

3. Dust in AGNs and AGN power indicator 

First of all, we develop a basis to characterize AGN activi- 
ties. We first quantify dust extinction in AGNs. We then study 
emission line properties using the extinction corrected lumi- 
nosities and examine indicators of the AGN power Finally, we 
revisit the issue of measuring SFRs of the host galaxies from 
emission lines. 

3.1. Extinction in AGNs 

AGNs can be obscured by dust around the nuclear regions 
and it is essential to correct for the dust extinction in order to 
quantify AGN activities. We measure emission line luminosi- 
ties due to AGNs in a two-step procedure. First, we subtract 
the star formation component, Lsf, from a measured emis- 
sion line luminosity. We then correct for the dust extinction 
to obtain emission line luminosities coming out of the central 
region. To be specific, in case of Ha, 

Lna = (iffQ.meas " LHa,sf) X exp(0.3 X 0.75 X Ty), (2) 

where Lna.meas, Lna.sf and Ty are the measured Ha lumi- 
nosity, 'predicted' Ha luminosity due to star formation, and 
the optical depth in the y-band from the spectral fitting. The 
factor of 0.3 to ry is due to the assumption that 30% of the 
dust extinction comes from the ambient interstellar medium 
(Chai-lot & Fall 2000). The term G.75Ty is the optical depth 
at the wavelength of Ha assuming the Cardelli et al. (1989) ex- 
tinction curve. We note that we tend to underestimate ry from 
the spectral fitting (see Fig. 1 of Paper-I) and the extinction due 
to the interstellar medium may not be fully corrected'. 

We apply the same equation to Lh^ and plot the Balmer 
decrement of the AGN emissions in Fig. 1. Here we adopt 

' We underestimate ry systematically by Ary ~ 1 for most objects as 
shown in Fig. 1 of Paper-I. If we empirically correct for the offset, we 
obtain a median AGN extinction of ry =0.7, which is smaller than that 
obtained below (ry = 1.2). 
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the dust-free B aimer decrement of 3.1 for AGNs (Halpern 
& Steiner 1983; Ferland & Netzer 1983; Gaskell & Ferland 
1984). We also have confirmed that there is no change in our 
main conclusions if we adopt a decrement of 2.86 (case 'B' re- 
combination for T = 10^ K and = 10^ cm ''; Osterbrock 
& Ferland 2006). The Oxygen-excess objects typically have 
Ty = — 2 extinction with the median being Ty = 1.2. We note 
that the median value does not change if we include lower sig- 
nificance detections of the lines (e.g., 3ct), although the scatter 
significantly increases. Dahari & De Robertis (1988) showed 
that Seyfert 2's typically show the Balmer decrement of 3 — 6, 
which is fully consistent with our estimates here. Kauffmann 
et al. (2003) quoted a median extinction of Ay ^ 3 mag. in 
AGNs using data from SDSS. Our extinction is less than that, 
but their large value is due to the assumption of the relatively 
flat extinction curve of Chariot & Fall (2000). We have con- 
firmed that we obtain a similar amount of extinction if we use 
the Chariot & Fall (2000) curve. It is interesting to note that the 
observed range of Balmer decrement is similar to that observed 
in Seyfert galaxies (Osterbrock 1977; Greene & Ho 2005). 

For the Balmer decrement measurements, we need good 
signal-to-noise ratios on the Balmer lines, but H/3 can often 
be noisy (and that is the reason why the BPT diagnostics is 
not applicable for a large number of galaxies). It is difficult 
to measure the Balmer decrement for individual galaxies and 
we therefore apply a constant extinction correction regardless 
of the properties of the AGNs and host galaxies. We note that 
we do not observe strong dependence of the balmer decrement 
on AGN power and thus the constant correction is a reasonable 
first order correction. The total optical depth of Oxygen-excess 
objects is 



TV,total = 0.3 X ry + 1.2. 



(3) 



The first term is the extinction due to the ambient interstellar 
medium, which is derived from the spectral fitting. The second 
term is the extinction around the nuclear region. To be specific, 
we replace 0.3 x Ty in Eq. 2 with Eq. 3 to obtain intrinsic line 
luminosities due to AGN. We use a subscript AGN when we 
show star formation subtracted, extinction corrected emission 
line luminosities of the Oxygen-excess objects. This is just for 
simplicity and we do not mean that all of the Oxygen-excess 
objects are real AGNs. 

3.2. A new indicator of AGN power and comparisons with 
X-ray power 

Having quantified extinction in AGNs, we are now ready to 
characterize AGN activity. In order to quantify AGN power, 
we ideally want to integrate all the photons emitted from the 
central engine over the whole wavelength range. But, this is 
observationally very challenging to achieve. A luminosity of 
the [Olll] line is often used as an indicator of AGN power and 
a clear correlation between the absorbed [Olll] luminosity and 
hard X-ray luminosity has been observed in a hard X-ray se- 
lected sample (Heckman et al. 2005). Recently, de Gasperin et 
al. (2011) also presented a tight correlation between absorbed 
hard X-ray luminosity and [oill] luminosity. The use of [Olll] 
as a proxy for AGN power is partly motivated by a commonly 
adopted assumption that contamination from star formation is 
relatively small in [oill]. But, [Olll] is a relatively high ion- 
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Fig. 2. Fraction of [Olll] agn luminosity to [OII]yi(jjv+[OIII] AGiV 
luminosity plotted against [OII].4(3jv+[oni]^Gjv- The solid contours 
are Oxygen-excess objects and their star formation fluxes are sub- 
tracted off and the extinction is corrected for using Eq. 3. The hor- 
izontal dashed line shows the equality of [011]agn and [omlylGiV 
luminosities. [Oil] agn<[OIU] AGN holds above the line, and 
[Oil] AGN > [Olll] AGN holds below the line. 

ization line. Is it a good indicator of AGN power for low- 
luminosity AGNs studied here? 

If we look at ionization state of the Oxygen-excess galax- 
ies in Fig. 2, we find that AGNs spans a wide range in ion- 
ization state. Heckman (1980) first introduced a class of low- 
ionization emission line galaxies (LINERs). The original def- 
inition of LINERs is [OII]>[OIII] and [OI]> 0.33[OIII]. 
Ho et al. (1993) showed that this definition is essentially 
the same as [OIII]/H/3 > 3. Due to the weak H/3 emis- 
sion in most of the Oxygen-excess galaxies, we further re- 
vise the definition to [OII]yiGjv>[OIII]^Gjv, which is equiv- 
alent to [OIII]/H/3 > 3 due to the tight correlation between 
[OII]/[OIII] and [OIII]/H/3 (Baldwin et al. 1981). Note 
that we use star formation subtracted, extinction corrected lu- 
minosities for the LINER/Seyfert classification. At low lu- 
minosities, AGNs are mostly LINERs. Seryfert-like galax- 
ies with [OII]yiGA'<[OIII]AGJV appear at high luminosi- 
ties. We find that 7% of the Oxygen-excess objects have 
[Oll]yiGJV<[Olll]AGJV- This fraction should not be taken as 
a general fraction of Seyferts among AGNs because we inten- 
tionally excluded QSO-like objects from the sample. 

[oill] is often used as an indicator of AGN power (Heckman 
et al. 2005), but given the relatively wide range of ionization 
states of the Oxygen-excess objects, we deem that the sum of 
[OII]h-[OIII] is a better indicator of the AGN power. In gen- 
eral, [on] and [Olll] lines balance each other in gaseous nebu- 
lae, [oill] is strong in high ionization states and [oil] is strong 
in low ionization states. A sum of the two lines is more ro- 
bust to variations in ionization levels than either one of them 
and is a better indicator of total energy emitted by Oxygen. 
Ha line would be another straightforward option. However, as 
shown in Paper-I, AGNs exhibit relatively weak Ha and the 
AGN component does not stand out. If we subtract the star 
formation component as we do for [Oil] and [oill], we ob- 
tain a negative Ha luminosity for a number of Oxygen-excess 
objects. This means that our SFR and ry estimates from the 



Tanaka 



5 



spectral fits are not sufficiently good to extract the relatively 
weak Ha due to AGN. On the other hand, the collisionally ex- 
cited lines are normally strong in AGNs. For these reasons, 
we choose [oil] + [oill] to identify AGNs and to characterize 
them. If one could improve the estimates of Ha flux due to star 
formation. Ha would be a better indicator because it is less 
sensitive to ionization, dust, and metallicity. 

We now test this new indicator of the AGN power by com- 
paring with hard X-ray luminosities, which is likely a very 
good indicator of the AGN power (Ho 2008). As described 
in Paper-I, we use the Chandra-SDSS matched catalog release 
1.1 (Evans et al. 2010). We require the sources to be point-like 
and located within 2 arcsec from the optical center to ensure 
that X-rays likely originate from the center of the galaxies. We 
compare the optical emission line luminosities with X-ray lu- 
minosities to verify the effectiveness of the new AGN power 
indicator We note that we remove X-ray sources that may 
be non-AGN origin from the comparisons (Fig. 11 of Paper- 
I show how we identify those sources). We further note that 
we do not correct for absorption in hard X-rays. The typical 
Galactic neutral hydrogen column density in our sample is low 
(a few times 10^'' cm^^; Dickey & Lockman 1990) and intrin- 
sic absorption is likely below lO^^cm^^ within the luminos- 
ity range we explore here (Mainieri et al. 2007). At this col- 
umn density, absorption in the hard band is small (Morrison & 
McCammon 1983). Our hard X-ray luminosity should there- 
fore be a reasonable estimate. 

We show comparisons between hard X-ray luminosity and 
optical emission line luminosities in Fig. 3. We show not only 
[oii]h-[oiii]agw, but [oii]agw and [oiiiI^gat as well. We 
find that all the lines show statistically significant correlations 
with X-rays. We perform the Spearman's rank test and find 
the positive correlation with a high significance (> 2a) to re- 
ject the null hypothesis in all the lines. Also, we perform the 
biweight fits to the data points and show the best-fitting lines 
in the plots. The dominant source of error in our AGN flux 
estimates is in our prediction of emission line fluxes due to 
star formation, which is 0.23 and 0.29 dex for [oil] and [oill], 
respectively. We adopt it as an uncertainty in the subtraction 
of the star formation component for each line. The biweight 
fits are performed taking these uncertainties into account. We 
find that [oiiiJ^gtv shows a stronger correlation with the X- 
ray power than [Oll]AGAf (Simpson 1998). Contrary to our 
expectations, [oil]H-[Olll]yiGJV shows somewhat poorer cor- 
relation than [oill]yiGJV alone, although the scatter decreases 
slightly. Apparently, the relatively poor correlation between 
{011\agn and X-ray luminosities weaken the correlation in 
[Oll]H-[Olll]^gjv We suspect that it might be due to insuf- 
ficient extinction correction. Due to our limited ability to mea- 
sure the extinction in individual galaxies, we apply the constant 
correction for dust due to AGN to all the galaxies. But, this is 
probably just a rough correction and a small over/under cor- 
rection introduces a larger error in [oil] than in [Olll] because 
[on] is at a shorter wavelength. 

Now, we face a difficult question of the choice of the AGN 
power indicator. We would like to use [oil]H-[oill] to keep 
consistency with the identification method. As shown in Fig. 
7 of Paper-I, the stacked spectra of the Oxygen-excess objects 
show stronger [oil] than [Olll] and [Oil] is a very powerful 



-1 -0.5 0.5 1 

Hardness ratio 

Fig. 4. [oni]AG]v/([OII]AG]V+[OIII]AGiv) plotted against X-ray 
hardness ratio between soft (0.5-1.2 keV) and hard (2-7keV) bands. 
The numbers show a Spearman's rank conelation coefficient and a null 
probability. 

feature to identify Oxygen-excess objects. On the other hand, 
we would like to use [oill] only as a better indicator of the 
AGN power at the cost of loosing consistency with the identifi- 
cation process somewhat. We have performed all the analyses 
presented in the paper using the two AGN power indicators 
and found that our conclusions remain totally the same regard- 
less of the choice of the indicator. Given this robustness of 
our conclusions, we choose [oil]H-[oill] as an indicator of the 
AGN power to keep the consistency with the AGN identifica- 
tion scheme. 

We note in passing that we find an interesting correlation be- 
tween "optical hardness" and X-ray hardness in Fig. 4. The 
optical line ratio plotted is essentially a simple function of 
[oil]/[oill], which is tightly correlated with [oill]/H/3 and is 
a good indicator of the ionization level (Baldwin et al. 1981). 
The Spearman's rank correlation coefficient suggests a posi- 
tive correlation between the optical and X-ray hardness with a 
statistically significant probability to reject the null hypothesis. 
This correlation can be understood as a product of two other 
correlations: one between Lx and hardness, and the other be- 
tween emission line luminosity (which is correlated with Lx) 
and ionization state. We find that most of the hard sources in 
Fig. 4 are powerful AGN with Lx > lO"^-^ erg s^^. The AGN 
power then correlates with the ionization state as we will show 
in Section 4. These two correlations result in the trend in Fig. 
4. 

3.3. Possible evidence against a dominant role of post-AGB 
photo-ionization 

Let us go back to the issue of objects with non-AGN photo- 
ionization. As we discussed in Paper-I and earlier in this pa- 
per, photo-ionization due to post-AGB stars can be important 
in LINERs. This comes from a few observational pieces of ev- 
idence. One is an energy budget argument. AGNs in LINERs 
tend to show a deficit in the amount of ionizing photons to 
fully explain the observed emission line luminosity (Maoz et 
al. 1998; Eracleous et al. 2010), requiring an additional source 
of ionization mechanism other than AGN. Another piece of 
evidence is that spatially extended emission is observed in 
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Fig. 3. Optical emission line luminosities plotted against hard X-ray luminosity. The plots show [OU]agn, [OIII] 4G'jv. and [on]+[OIII]AGiV from 
left to right. We remove soft low-Lx sources, some of which may not be AGNs, from the plots. The numbers in each plot show a Spearman's rank 
coiTelation coefficient and a null probability. The solid lines ai'e biweight fits to the data points and the dispersion in the optical luminosity ai'ound them 
are also shown in the plots. 



LINERs using 2D spectrograph (Sarzi et al. 2010). Recently, 
Yan & Blanton (2011) also detected extended emission by a 
clever analysis of extended emission using varying redshift 
slices. We suspect that they suffer from contamination of star 
forming galaxies because a cut in D4000 is not enough to elim- 
inate them from their sample as we mention below. Extended 
emission does not necessarily support the post-AGB origin of 
the emission (Yan & Blanton 2011), but spatially extended 
post-AGB stars remain a natural explanation of it. It is not 
clear whether post-AGB stars alone can produce enough ion- 
izing photons to fully explain the observed emission (e.g., Yan 
& Blanton 201 1 reports on the deficit of ionizing photons even 
if they assume that 100% of the ionizing photons are absorbed 
to produce emission line, while Sarzi et al. 2010 claims that 
they have enough ionizing photons). But, there is no doubt 
that emission due to post-AGB stars contaminate the emission 
line luminosities observed in the Oxygen-excess objects. The 
question is how much post-AGB stars contribute to the overall 
emission. To address this question, we focus on the correla- 
tion between optical emission line luminosities and hard X-ray 
luminosities in Fig. 3. 

Post-AGB stars are not luminous in hard X-rays and the ob- 
served hard X-ray luminosity is likely primarily due to AGNs. 
One may worry that low-mass X-ray binaries may contribute 
to the hard X-ray luminosity. But, we show in Fig. 5 that 
the hard X-ray luminosity does not correlate with host galaxy 
stellar mass. The Spearman's rank correlation coefficient is 
very small and the test does not reject the null hypothesis. This 
rather weak (or no) dependence of X-ray luminosity and host 
galaxy mass has also been noted by other authors (Mullaney et 
al. 2012; Aird et al. 201 1). This strongly rules out the low-mass 
X-ray binary origin of the hard X-ray (Kim & Fabbiano 2004). 
The observed hard X-ray luminosity is primarily due to AGNs. 

On the other hand, optical emission can be contaminated by 
post-AGB photo-ionization. The correlation between optical 
emission line luminosity and hard X-ray luminosity can there- 
fore be used to gauge the contamination from post-AGB stars. 
The observed clear correlation between them shown in Fig. 3 
suggests that the post-AGB contamination is not significant. 
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Fig. 5. Hard X-ray luminosity plotted against stellar mass. Objects 
that are likely contaminated by non-AGN sources (those in the dashed 
box in Fig. 11 of Paper-I) are removed. We also show the Spearman's 
rank congelation coefficient and null probability. 

If the contamination was significant, we would not have ob- 
served such a clear correlation between hard X-ray and, e.g., 
[oil]-i-[oill] luminosities. We shall note that the X-ray sources 
are drawn from serendipitous Chandra observations and they 
have typical properties of the Oxygen-excess objects, i.e., they 
are massive galaxies with > lO^^M© in stellar mass and are red 
quiescent galaxies with SFRs below O.IMq yr~^. Therefore, 
these X-ray sources should be representative of the Oxygen- 
excess objects. 

The X-ray sample is a small subset of the entire Oxygen- 
excess objects. To further quantify the emission due post-AGB 
stars, we define a subsample of the Oxygen-excess objects lo- 
cated at 0.055 < z < 0.060 and with zero SFRs from the spec- 
tral fits. The narrow redshift slice is adopted to minimize the 
variation of our sensitivity to emission line luminosities within 
the redshift slice. The redshift range is arbitrarily chosen, but 
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Fig. 6. [OII]+[oni] luminosity plotted against stellar mass. The 
Oxygen-excess objects plotted here are located at 0.055 < z < 0.060 
and are not forming stars. Extinction defined in Eq. 3 is not applied and 
thus the luminosity is a 'raw' luminosity. Stellar mass is not corrected 
for the aperture loss. The solid line is a biweight fit to the data with a 
power-law slope of 0.26. The dashed curve shows the best fit model to 
the data (see text for details). 

our conclusion below does not strongly depend on a particular 
choice of the redshift range. The zero SFR is to ensure that 
the host galaxies are not forming stars so that we do not have 
to correct for emission line luminosities due to star formation. 
The observed emission is therefore due to AGN and to post- 
AGB star photo-ionization (if present) in this subsample. 

We show [oil]-i-[oill] luminosity against stellar mass of this 
subsample in Fig. 6. Here, we do not apply any dust correc- 
tion or aperture correction, i.e., 'raw' emission line luminosity 
within the fiber aperture is plotted against stellar mass within 
the fiber. The emission due to post-AGB stars should corre- 
late strongly with stellar mass contained within the same area 
covered by the fibers. On the other hand, the AGN emission 
is unlikely to be strongly correlated with mass. Therefore, the 
correlation in Fig. 6 is a rough measure of the contribution of 
the post-AGB stars. 

The Spearman's rank test suggests that there is a correlation 
between the emission line luminosity and stellar mass. The 
null hypothesis is rejected at a significant level. This suggests 
that the post-AGB stars do contribute to the observed line lumi- 
nosities. However, the observed correlation is relatively weak 
- the biweight fit gives a log-linear slope of 0.26. That is, if 
stellar mass increases by an order of magnitude, the emission 
line luminosity increases by less than a factor of 2. The abun- 
dance of post-AGB stars should depend on star formation histo- 
ries of galaxies, which then depend on galaxy mass. However, 
we are looking only at massive quiescent galaxies-, and it is 
very unlikely that their star formation histories vary so widely 
that the stellar mass dependence of the post-AGB abundance is 
smoothed out. The most naive interpretation of Fig. 6 would 
be that post-AGB stars do contribute to the observed line emis- 

^ Stellar mass plotted in Fig. 6 is not corrected for the aperture loss and 
is significantly smaller than a typical stellar mass of the Oxygen-excess 
objects discussed later in this paper, which is corrected for the aperture 
loss. 



Fig. 7. SFR plotted against Dn,4000- The contours are all the objects 
in our parent sample (i.e.. Main galaxies at 0.02 < z < 0.10). The 
points and enor bars show the median and quailile of the SFR distribu- 
tion. 

sion, but their contribution is minor 

Let us make an attempt to quantify the contribution of post- 
AGB stars using a very simple model. If we assume that emis- 
sion due to post-AGB stars linearly increase with stellar mass 
and emission due to AGN does not depend on stellar mass, we 
can write 

Lobs = Lag N + Ct stellar, (4) 

where Lobs, Lagn, ot are observed luminosity, luminosity due 
to AGN, and luminosity due to post-AGB photo-ionization per 
unit stellar mass. Let us modify the second term to define post- 
AGB emission relative to Lagn at M stellar = 1O^°M0. 

T T / 1 , ^stellar \ 

Lobs^LAGN(l + P^J^^). (5) 

If /3 = 1, AGN and post-AGB equally contribute to the observed 
emission at Msteiiar = 1O^°M0. We allow Lagn and f3 to vary 
and fit the observed data in Fig. 6 using the statistics. We 
find logLAGN = 39.707 ± 0.005 and /3 = 0.305 ± 0.011 give 
the best fit. The errors are derived by taking Ax^ = 1 and are 
very small, but we should note that we ignore intrinsic scatter 
in Lagn- The best-fit model is shown as the dashed curve in 
the figure. The fit suggests that 23% of the observed emission 
is due to post-AGB stars at Msteiiar = lO^^M©. Because of the 
assumptions we made, this is nothing more than a rough guess. 
But, it seems that post-AGB stars do not play a dominant role. 

This result may not appear consistent with those from the 
literature (e.g., Sai'zi et al. 2010; Yan & Blanton 201 1). We do 
not try to resolve the issue as the SDSS data is unlikely a right 
data set to address it as we mentioned in Paper-I, but we shall 
recall that our concern is in a very low luminosity regime. We 
have to be sure that we do not suffer from low-level star forma- 
tion that could happen in early-type galaxies. To illustrate its 
importance, we show SFR against £'„,4ooo in Fig- 7. DnAooo 
is a sensitive index to age of the stellar population and is used 
in earlier work to select quiescent galaxies. However, it is hard 
to construct a pure sample of quiescent galaxies with this index 
alone. Even if galaxies have a large L'„,4oooi they may well be 
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undergoing low-level star formation as shown in the figure. A 
SFR of O.OIMq yr^i translates into Lua = 2 x lO^^erg s^^ 
(assuming no extinction by dust), which is a level of emission 
line luminosity that previous work reported. Such star forma- 
tion is likely spatially extended. A careful removal of low-level 
star formation would be necessarily to put a better constraint on 
the role of post-AGB stars. 

To summarize, emission line luminosities due to post-AGB 
photo-ionization seems to be a minor component of the over- 
all luminosities. However, other authors claim a substantial 
role of post-AGB stars (Sarzi et al. 2010; Cid Fernandes et al. 
2011; Yan & Blanton 2011). This issue remains unresolved and 
further effort is needed to settle it. One difficulty that hampers 
a clear conclusion is our limited understanding of the last stage 
of the stellar evolution. Another difficulty is the lack of suitable 
observational data to address the issue. Dedicated observations 
would be essential. For these difficulties, we cannot try to pur- 
sue the issue further in this paper, but we should bear in mind 
the potentially non-negligible contamination from post-AGB 
stars. Given the unclear (but likely minor) role of post-AGB 
stars, we do not correct for their contamination to the observed 
emission line luminosities throughout the paper. We assume 
that the observed emission is entirely due to AGN and SF. 

3.4. Comparisons with radio power 

We turn our attention to radio sources. We cross-matched 
our objects with FIRST (Becker et al. 1995; White et al. 1997) 
within 1 arcsec apertures as described in Paper-I to compare 
optical power with radio power We miss extended radio 
lobes/jets, but extended radio sources constitute only 10% of 
the radio galaxies and hence should not strongly affect our re- 
sults here (Lin et al. 2010). Also, extended radio lobes may 
not be a good measure of instantaneous AGN activity. We use 
fintegrated from the FIRST Catalog. Our results remain the 
same if we use f peak- 

We plot the optical power against radio power in Fig. 8. 
Note that we remove radio sources that are dominated by star 
formation from the comparisons (Fig. 13 of Paper-I shows how 
we identify star formation dominated sources). We observe a 
weak correlation between optical and radio and the Spearman's 
rank test rejects the null hypothesis. The weak correlation is in 
contrast to X-rays, which show stronger correlations. The loose 
radio-optical correlation may appear inconsistent with previous 
studies (e.g., Saunders et al. 1989; Baum & Heckman 1989; 
RawHngs & Saunders 1991; Hes et al. 1993; Ho 1999). These 
authors observed a relatively tight correlation between optical 
emission and radio power over a range of AGN types. 

We do not consider our finding inconsistent with earlier find- 
ings for two reasons. First, our radio galaxy sample has more 
than an order of magnitude lower radio power than most of the 
previous studies. Majority of the Oxygen-excess objects are 
low-luminosity AGNs, while the previous studies focused on 
powerful radio galaxies. We are therefore not looking at the 
same population. Second, we exclude sources that are strongly 
contaminated by underlying star formation activities, while the 
previous studies did not. We find that the correlations slightly 
tighten up if we include those sources because we include the 
strong SFR-radio power correlation in the statistics. Our re- 
sult may still not be consistent with observations by Ho (1999), 



who studied nearby low-luminosity AGNs in early-type galax- 
ies, in which the underlying star formation is probably not very 
strong. However, their early-type galaxy sample shows a large 
scatter comparable to ours (see their Fig. 2). For these rea- 
sons, we do not consider that our results are in strong conflict 
with previous studies. Recently, de Gasperin et al. (2011) ob- 
served no clear correlation between radio and X-ray powers. 
This probably suggests that there is no strong correlations be- 
tween radio and optical powers in consistent with our result 
here. 

The positive correlation shows the black hole activities and 
radio activities are weakly correlated in the Oxygen-excess ob- 
jects under study. We note that several objects are detected both 
in X-rays and radio. We compare the X-ray and radio pow- 
ers and find no significant correlations. Our statistics is very 
poor, but the observed lack of correlation is consistent with de 
Gasperin et al. (201 1). The optical and X-ray powers are likely 
a good indicator of on-going AGN activities, while the radio 
power is likely not. 

3.5. Is [oil] an indicator ofSFRs of AGN host galaxies? 

The [on] emission is sometimes used as a proxy for host 
galaxy SFRs of high ionization AGNs such as Seyferts because 
this line is not very strong in those AGNs (e.g., Ho 2005). It 
may be reasonable to assume that [Oil] is due to star formation 
in Seyfert galaxies, but is it still good for LINERs? The newly 
developed method allows us to statistically separate star forma- 
tion and AGN components of a given emission line. We take 
this opportunity to quantify how well [oil] is correlated with 
SFRs over a range of AGN power and ionization states. 

Fig. 9 shows the fractional contribution of [oil] due to 
star formation to the overall [oil] luminosity as a function 
of SFRs of the host galaxies (left panel), AGN power (mid- 
dle panel), and ionization level (right panel). We do not 
have a lot of Oxygen-excess galaxies dominated by star for- 
mation (i.e., those above the dashed line), but this is sim- 
ply a selection bias that we cannot find weak AGNs in ac- 
tively star forming galaxies because our method uses a con- 
trast between Lsf and Lagn- If [oil] was a perfect indica- 
tor of star formation, we would expect that galaxies line up at 

L[oii],sf/{L[oii],sf + L[oii],agn) = 1- 

Starting with the left panel, [oil] due to AGN completely 
dominates in galaxies with low SFRs, while the star forma- 
tion component comes in at moderate-high SFRs. This is a 
totally expected behavior. The observed [oil] does not give 
any sensible SFRs in quiescent galaxies. Even in actively star 
forming galaxies {SFRspecfu ^ 1 yr~^), [Oil] due to star 
formation does not dominate, although we are limited by the 
selection bias there. 

One might expect that the AGN component becomes neg- 
ligible in weak AGNs and [oil] might be a good indicator of 
SFRs. But, the middle plot in Fig. 9 shows that this is not the 
case. As we will discuss in the next section, there is a positive 
correlation between the AGN power and underlying star for- 
mation activities. Due to this correlation, there is no strong de- 
pendence between the dominance of SF/AGN component and 
AGN power. There is always a fair chance of suffering from 
AGN contribution in the observed [Oil] luminosity. 

Finally, the right panel shows the dominance of the star for- 
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Fig. 8. As in Fig. 3, but liere we compare optical emis.sion line luminosities with radio power FIRST sources dominated by star formation are removed 
from the plot. The median en'or is shown on the bottom-right corner. The integrated FIRST fluxes are measured by fitting elliptical Gaussian to the 
sources, which makes it difficult to evaluate the flux uncertainties (White et al. 1997). We do not quote the median uncertainty in the radio power. 




Fig. 9. Left: Fractional contribution of [Oil] due to star formation to the overall [Oil] luminosity plotted against SFR. The dashed line shows the equality 
of SFR and AGN components. Star formation dominates above the line and AGN dominates below the line. The points with en'or bars show the median 
and quartiles of the distribution. Middle: Same as the left panel, but here the horizontal axis is AGN power. Right: Same as the left panel, but the 
horizontal axis is ionization state. [OII]<[OIII] holds at [OIII]/([OII]+[OIII])> 0.5, and [OII]>[OIII] holds at [OIII]/([OII]+[OIII])< 0.5. 



mation component as a function of ionization level. IVIost of 
the LINERs {[011\agn>[0111\agn) are dominated by AGN 
emission and [oil] due to star formation dominates only in the 
very high ionization galaxies. For the vast majority of low- 
luminosity AGNs, [on] is unlikely a good indicator of star for- 
mation. 

To be further quantitative, let us ask a simple question — 
if we assume that [oil] is entirely due to star formation, how 
well can we reproduce the SFRs of the hosts? Left panel in 
Fig. 10 gives the answer We simply translate an observed 
[on] luminosity into SFR using the formula given in Kennicutt 
(1998) without any subtraction of the AGN component or any 
dust correction. As can be seen, the SFRs from [oil] do not 
correlate with true SFRs very well. In particular, one grossly 
overestimates SFRs in quiescent galaxies. In those galaxies, 
star formation is very weak and the observed [oil] is nearly 
entirely due to AGN (Fig. 9 left panel). The assumption of 
[on] due to star formation completely fails here. Another point 
to mention is that the systematic offset seen at ^ IM0 yr^^ is 
mostly due to the ignorance of dust. If we correct for ry = 1 
extinction, which is typical in star forming regions (Hopkins et 
al. 2003), the systematic offset is largely gone. However, the 
dust correction increases the discrepancy in low SFR galaxies 
because it shift all the points upwards. 

As shown in Paper-I, the collisionally excited [Oil] and 
[Olll] suffer (or benefit when we identify AGNs using these 



lines) from AGN contribution, while Ha is less affected. This 
is consistent with what the BPT diagram tells us — AGNs show 
stronger collisionally excited lines than the Balmer lines and 
that is why galaxies form a distinct sequence from the star for- 
mation sequence to the top-right corner of the diagram, not 
to the bottom-left corner. Due to the relatively weak AGN 
contribution to Ha, SFRs from Ha indeed better agree with 
those from the spectral fits as shown in the right panel of Fig. 
10. Therefore, Ha is a better indicator of SFRs in AGNs than 
[on]. However, Ha still largely overestimates SFRs in quies- 
cent galaxies by more than an order of magnitude. 

To sum up, one should be cautious about using [oil] as an 
indicator of SFRs in AGNs. It can give SFRs that are offset 
by more than an order of magnitude. In particular, we sug- 
gest that one should not use [oil] when a host galaxy has red 
colors. The galaxy may be a dusty galaxy, but it may well be 
passive galaxies, for which [oil] is entirely due to AGN. As 
a result, one wildly overestimates SFRs. Also, dust extinction 
needs to be accounted for One should use continuum infor- 
mation whenever available because the spectral fitting is likely 
to deliver better SFRs due to less AGN contamination to the 
continuum spectrum (Schmitt et al. 1999). It may often be the 
case when the continuum spectrum is not available with suffi- 
cient signal-to-noise ratios especially at high redshifts. Recent 
availability of multi-wavelength photometry would give a way 
around it. SED fits of multi-wavelength broad-band photome- 
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Fig. 10. SFRs from [Oil] (left) and Hq (right) plotted against SFRs from the spectral fits. We only correct for the fiber aperture loss and do not apply 
any further connection to [Oil] and Ho (i.e., no subtraction of AGN component and no correction for dust) and simply convert the observed luminosity 
into SFRs. The points with en'or bars show the median and quartiles of the distribution. 



try seem to give reasonable SFR estimates (e.g., Santini et al. 
2009). If one is careful enough to remove photometry that can 
be affected by AGN component (UV continuum or IR thermal 
emission) from the SED fits, it would provide better SFR esti- 
mates than [on] and Ha. We shall emphasize that the key is to 
use stellar light only. It would be challenging to estimate SFRs 
of QSOs, whose continuum is dominated by AGN emission, 
not by stars. If QSO continuum can be subtracted by assuming 
a power-law form, one can fit the spectrum of the remaining 
stellar component. An accuracy of SFRs from such a fit re- 
mains to be quantified. 

4. Nature of the Oxygen-excess objects: AGN activities 
and the host galaxy properties 

Following the development of the basis to characterize AGN 
activities, we now move on to discuss the nature of Oxygen- 
excess objects with an emphasis on the relationship between 
the AGN and host galaxy properties. We start with the fraction 
of AGNs as a function of host galaxy properties, followed by a 
detailed look at a partition between LlNERs and Seyferts. We 
then show scatter plots of the host galaxy properties and black 
hole masses. Finally, we study the AGN activity and compare 
the black hole growth rate and the host galaxy growth rate. 

4.1. Dependence of AGN fraction on Stellarmass, Color, SFR 
and Morphology 

Essentially all galaxy properties depend on mass of galaxies 
and we start with stellar mass dependence of the AGN frac- 
tion in Fig. 11. The AGN fraction is a very strong function of 
stellar mass of host galaxies. Most AGNs that we have identi- 
fied reside in galaxies with Msteiiar,apercorr > 1O^°M0. In the 

most massive galaxies with Mstellar,apercorr > IO^'^^Mq, 60% 

of the galaxies harbor AGNs. We note that this is just a lower 
limit. A large fraction of galaxies in our sample do not exhibit 
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Fig. 11. Fraction of AGNs plotted against stellar mass. The filled and 
open circles are Oxygen-excess objects and AGNs from the BPT dia- 
gram, respectively. The solid line shows the fractional distribution of 
all the galaxies in our sample. 

Strong enough Oxygen lines to apply our method (i.e., OnBn 
defined in Paper-I) and the stacked spectrum of those galaxies 
exhibits LINER-like line ratios. Therefore, many of the un- 
classified galaxies may well be AGNs (see also discussions on 
the non-AGN photo-ionization in Paper-I and earlier sections 
of this paper). In contrast to the Oxygen-excess method, the 
BPT diagnostics gives an AGN fraction of only 20% in mas- 
sive galaxies. This demonstrates the sensitivity of our method 
to identify AGNs. 

Heckman et al. (2004) found that the fraction of active black 
holes (i.e., AGNs) is highest in intermediate mass black holes 
and it decreases in the most massive black holes based on BPT 
AGNs. We have replaced the stellar mass with black hole mass 
(see below for the derivation of the black hole mass) in Fig. 1 1 
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and found that a fraction of Oxygen-excess objects continues 
to increase towards the high mass end. This difference from 
Heckman et al. (2004) is due to the difference in the redshift 
ranges explored. Heckman et al. (2004) did not apply any red- 
shift selection to the main galaxy sample (Strauss et al. 2002). 
Due to the flux limited nature of the sample, an average red- 
shift of galaxies is higher for more massive galaxies and there- 
fore an average sensitivity to emission line luminosity is lower 
This resulted in a reduced sensitivity to identify AGNs in very 
massive galaxies and caused a decline in the fraction of active 
black holes at the high mass end. We have confirmed that we 
reproduce the Heckman et al. (2004)'s trend by changing the 
redshift range to < z < 0.3 and redoing the whole analysis. 

We further discuss color, SFR and morphology dependence 
of the AGN fractions. 

Color (Fig. 12): We compute rest-frame u — r color from 
the model magnitudes (Stoughton et al. 2002) using the k- 
correction code by Blanton & Roweis (2007) and plot the 
AGN fraction as functions of stellar mass and rest-frame u — r 
color. We find that Oxygen-excess objects prefer red galaxies 
(w — r ^ 2) at all masses. As expected from the stellar mass 
dependence in Fig. 11, the fraction increases from low to high 
mass galaxies. The most massive galaxies show the highest 
AGN fraction, reaching to nearly 70%. The fraction decreases 
towards bluer colors, but this is at least partly driven by the se- 
lection bias that we cannot identify weak AGNs in actively star 
forming galaxies. 

In contrast to the clear preference of Oxygen-excess objects 
to red colors, a large fraction of BPT AGNs reside in the green 
valley around m — r 2 in consistent with Kauffmann et al. 
(2003) and Schawinski et al. (2010). This is due to the fact 
that BPT is not very sensitive to AGNs in passive galaxies. As 
shown later, most AGNs in passive galaxies are low-luminosity 
AGNs. They often show too weak H/3 to apply the BPT method 
(see the stacked spectrum of O-nBn in Fig. 7 of Paper-1). As a 
result, the BPT misses a significant fraction of low-luminosity 
AGNs on the red sequence, and that causes the large difference 
in the AGN fraction in red galaxies between Oxygen-excess 
and BPT. In other words, BPT AGNs in the green valley 
are relatively strong AGNs. We note that objects detected 
in radio or X-ray do not show a clear excess in the green valley. 

SFR (Fig. 13): In addition to the color, SFR is also one of 
the most important properties of galaxies. We quantify the 
AGN fraction as a function of SFRs and stellar mass of the host 
galaxies. Low-mass galaxies are predominantly star forming, 
but massive galaxies exhibit bimodal distribution of SFRs. The 
BPT AGN fraction shown in the right panel is constantly low 
at SFR ^ O.IM0 yr~^, but it shows a rapid increase at higher 
SFRs. On the other hand, the Oxygen-excess objects show 
a contrasting behavior. We interpret it by a combination of 
two effects. At low SFRs, the Oxygen-excess method is fairly 
sensitive and it gives a much higher AGN fraction than BPT. 
Because the method is based on identifying excess Oxygen lu- 
minosity, it becomes less efficient as underlying star formation 
becomes more active. This causes the drop of the AGN fraction 
around IM0 yr~^, where significant star formation component 
kicks in (see Fig. 9). This plot shows the sensitivity function of 



the Oxygen-excess method (we will further quantify it below). 
Its sensitivity is similar to BPT when SFRs are high, but it is 
much more sensitive at low SFRs. The Oxygen-excess method 
is suited to study AGNs in quiescent galaxies. 

The fraction of Oxygen-excess objects seems to more 
strongly depend on stellar mass than on SFR. The fraction of 
Oxygen-excess objects with SFR ~ 0.1 yr^^ drops from 
70% at high mass to nearly 0% at low mass. Mass is likely the 
primary parameter to trigger the AGN activity. 

Morphology (Fig. 14): As an indicator of morphology, 
we use the inverse concentration index (Shimasaku et al. 
2001; Strateva et al. 2001) defined by i?5o/i?9o, where i?5o 
and i?9o are the radii which enclose 50% and 90% of the total 
light (Petrosian flux; Petrosian 1976; Stoughton et al. 2002) in 
the z-band to minimize effects of on-going star formation. The 
inverse concentration index is of course affected by the seeing. 
The median seeing in the z-band is 1.3 arcsec. We have con- 
firmed that the trend in the plot does not change whether we 
use only objects observed under good seeing or not. 

We quantify the AGN fraction as functions of the inverse 
concentration index and stellar mass in Fig. 14. Oxygen- 
excess objects tend to be massive early-type galaxies as 
expected from the color-stellar mass dependence discussed 
above. The Oxygen-excess fraction shows a mild dependence 
on morphology, but again, stellar mass seems to be the primary 
parameter to trigger AGN activity. It is interesting to note 
that even if galaxies exhibit early-type morphology, they do 
not show a hint of AGN activity at low mass. In contrast 
to Oxygen-excess objects, BPT AGNs prefer intermediate 
morphology as they are strongly biased towards green valley 
galaxies due to the poorly sampled low-luminosity AGNs. 

Overall, Oxygen-excess objects tend to be massive red 
galaxies. The fraction of Oxygen-excess objects decreases to- 
wards blue, star-forming, late-type galaxies at a given stellar 
mass, but this is at least partly driven by the selection bias that 
we miss weak AGNs in actively star forming galaxies. On the 
other hand, the distribution of BPT AGNs is strongly skewed 
towards galaxies in the green valley. This is driven by the bias 
that BPT diagnostics is not very sensitive to low-luminosity 
AGNs on the red sequence. The figures presented in this sec- 
tion seem to suggest that BPT diagnostics does not provide an 
unbiased sample of AGNs. BPT AGNs are strongly biased. 
Oxygen-excess objects are also biased, but the selection bias is 
relatively straightforward to quantify as shown below because 
it is simply based on a contrast between AGN power and under- 
lying star formation (it can be complicated when the featureless 
continuum emission from AGN significantly contributes to the 
overall spectra). We will carefully account for such biases in 
the following discussions. 

4.2. Dependence ofLINERs and Seyferts on host galaxies 

Let us briefly mention the partition between LINERs 
and Seyferts among the Oxygen-excess objects as a func- 
tion of host galaxy properties. We define LINERs as 
[OIIIagw >[0111\agn as discussed in Section 3.2. We re- 
mind the readers again that the subscript AGN means that 
the luminosity due to star formation is subtracted and the ex- 
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Fig. 13. Left: Fraction of Oxygen-excess objects as functions of SFR and stellar mass. The contours are the same as in Fig. 12. Galaxies with 
SFRs< IO-^Mq yr-l are summed up at SFR= IQ-^M© yv"^. Right: As in the left panel, but for BPT AGNs. 
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Fig. 15. Fractions of Seyferts and LINERs as a function of stellar 
mass. The filled and open circles show LINERs and Seyferts, re- 
spectively. The solid line shows the normalized distribution of all the 
Oxygen-excess galaxies. 

tinction correction described in Section 3.1 has been applied. 
Kewley et al. (2006) and Schawinski et al. (2007) studied the 
LINER/Seyfert partition using BPT AGNs, but we revisit the 
subject with pure AGN luminosities. 

As in the last subsection, we start with stellar mass depen- 
dence shown in Fig. 15. It is striking that most Oxygen-excess 
objects are LINERs and the LlNER/Seyfert partition shows 
only very weak dependence on stellar mass. This is in con- 
trast to the strong stellar mass dependence of the AGN frac- 
tion discussed above. AGNs prefer massive galaxies, but once 
they become active, their activity (LINER vs Seyfert) does not 
strongly depend on mass. We will further discuss this 'mass 
switch' to activate AGNs in the next section. 

Fig. 16 shows the LINER/Seyfert partition as functions of 
stellar mass, rest-frame u — r color, SFR, morphology and 
AGN power There is a weak trend that Seyferts increase 
in late-type, blue, star forming galaxies. Massive red, quies- 
cent galaxies do not show Seyfert-like activities. This is qual- 
itatively consistent with Kewley et al. (2006) who observed 
that typical age of stellar populations in Seyferts is younger 
than that of LINERs. It is also in agreement with Schawinski 
et al. (2007) that LINERs tend to be in quiescent galaxies. 
The LINER/Seyfert partition is apparently not strongly de- 
pendent on host galaxy properties, but we find that the par- 
tition is a relatively strong function of the AGN power At 
L[oii]+[oiii]agn < 10"^^ erg s-\ the Oxygen-excess objects 
are mostly LINERs. The Seyfert fraction rapidly increases 
towards high luminosity and Seyferts and LINERs almost 
equally exist in the most powerful objects studied here. We 
note that, within the redshift range explored (0.02 < z < 0.1), 
we are complete down to a luminosity of L[oii]+[oiii]agn ^ 
3 X 10"'" erg for quiescent galaxies. This number is de- 
rived from the lower boundary of galaxy distribution on a 
L[oii]+[oiii]agn vs redshift plot. 

To summarize, the partition is primarily a function of AGN 
power and its dependence on host galaxy properties is rel- 
atively weak. Seyferts seem to slightly prefer star forming 
galaxies, but this correlation is likely due to a correlation be- 



tween SFR and AGN power presented below. 

4.3. Correlations between AGN activity and host galaxy 
properties 

To further discuss the AGN activity, we introduce black hole 
mass. There are a few ways to estimate black hole mass. Here 
we take an empirical way and use the stellar velocity disper- 
sion. Bulge stellar velocity dispersion is observed to exhibit 
a tight correlation with black hole masses (Gebhardt et al. 
2000; Ferrarese & Merritt 2000; Tremaine et al. 2002; Giiltekin 
et al. 2009). We use the velocity dispersion estimates from 
the MPA/JHU catalog (Tremonti et al. 2004; Kauffmann et al. 
2004) and derive the black hole masses using the empirical re- 
lation from Tremaine et al. (2002): 

log(MB///Afo) = 8.13 + 4.021og(a/200), (6) 

where a is the stellar velocity dispersion in units of 
km s~^. This calibration is consistent with a recent calibration 
(Giiltekin et al. 2009). We should emphasize that our black 
hole mass estimates are nothing more than a rough guess and 
are only of a statistical value as warned by previous studies 
(e.g., Kauffmann et al. 2007). The SDSS measures the velocity 
dispersion within 3 arcsec fibers, and hence it is measured over 
different physical sizes for galaxies at different redshifts. A ve- 
locity component of disk will be a significant source of error in 
our black hole mass estimates. There is no straightforward way 
to correct for such an effect and dedicated observations would 
be needed for more accurate black hole mass estimates. Note 
that we remove objects with cr < 70km from the catalog 
because of the limited SDSS spectral resolution of i? ~ 2000. 

We show in Fig. 17 correlations between the black hole 
mass and stellar mass, SFRs, morphology, and AGN power. 
Although not all the plots are entirely new and interesting, we 
show them as a representative set of correlations between the 
primary parameters. We make the following points from these 
plots. 

Stellar mass: Black hole mass increases with increasing 
stellar mass. This correlation is expected given that more 
massive galaxies show higher stellar velocity dispersions. We 
obtain a Spearman's rank correlation coefficient of 0.58 with a 
null probability of <C 10'^. Most Oxygen-excess objects have 
Mbh between 10^ and 10^ M©. This is similar to the black 
hole mass distribution of LINERs and Seyferts in the Palomar 
survey (Ho 2009). 

SFR: Given the clear positive correlation between stellar mass 
and black hole mass, this plot is basically equivalent to a SFR 
vs stellar mass plot. Galaxies form two well-known sequences 
- star forming and quiescent. The bimodality is particularly 
clear in massive galaxies. 

Morphology: Most of the Oxygen-excess objects exhibit 
early-type morphology with the inverse concentration index 
below 0.4 (Shimasaku et al. 2001; Strateva et al. 2001). 
The index gradually increases towards low black hole mass, 
but we are not sure if this is a real trend. Lower mass 
objects appear smaller on the sky and thus they are more 
strongly affected by the seeing effect, which results in in- 
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Fig. 16. Fraction of Seyfert galaxies ([Oil] ag N <lOUl] AG n) among all the Oxygen-excess objects. The panels show the fraction as functions of stellar 
mass and rest-frame u — r color (top-left), SFR (top-right), inverse concentration index (bottom-left), and AGN power (bottom-light). We show our 
sensitivity limit of AGN power in the bottom-right panel as the dashed line. The contours are as in Fig. 12. 



creased inverse concentration index. But, in any case, it is fair 
to say that most Oxygen-excess objects are early-type galaxies. 

AGN power: There seems to be a weak positive correlation 
between the AGN power and black hole mass, although the 
scatter is very large. A Spearman's correlation coefficient 
is 0.10, but the null hypothesis is rejected with a <C 10^^ 
probability. Such a weak correlation is surprising given the 
strong stellar mass dependence of the AGN fraction discussed 
in Section 4.1. We recall that we did not observe a clear 
correlation between hard X-ray luminosity and stellar mass in 
Fig. 5. A similar trend has been observed by other authors 
(Aird et al. 2011; Mullaney et al. 2012). We will elaborate on 
the trend in Section 5. 

4.4. The co-evolution of black holes and host galaxies 

We have quantified the basic correlations between the pri- 
mary parameters and we now characterize relationships be- 
tween activities and growth rates of AGNs and those of host 
galaxies. Fig. 18 shows the correlation between AGN power 



and host galaxy SFR. This plot therefore compares the activ- 
ities of the host galaxies in terms of star formation (Msteiiar) 
and the accretion rate of the central black hole (Mbh)- We 
recall that we have subtracted the star formation component 
from the [Oll]H-[Olll] luminosity and corrected for extinction, 
thereby representing intrinsic AGN emission. The slanted 
dashed lines are our selection bias mentioned several times by 
now. Because we require a significant Lagn with respect to 
LsF in Eq. 1, we miss weak AGNs in actively star forming 
galaxies. The dashed lines are the limits where Lagn can be 
detected at 1.5cr (> l.ha is the definition of Oxygen-excess) at 
a given Ls f for galaxies with typical extinction. Unlike BPT, 
the selection bias of the Oxygen-excess method is understood 
fairly well. 

The AGN power seems to increase with increasing SFRs. 
But, the trend is driven by the incompleteness at the bottom- 
right corner of the plots due to the selection bias. The plot 
does not give us an average AGN power as a function of 
SFR, and we can only discuss the upper envelope of the 
AGN power distribution. AGNs in low-SFR galaxies are typi- 
cally of low-luminosity and few powerful AGNs (several times 
10^^ erg s"^, say) are located in those galaxies. The upper 
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envelope of the AGN power distribution extends to higher lu- 
minosity at higher SFRs and powerful AGNs are often hosted 
by actively star forming galaxies. This suggests that galaxies 
and the central black holes can become active simultaneously. 
This is a surprising phenomenon because galaxy star formation 
and the central AGN activity occur on a very different physical 
scales (10 kpc vs. 1 pc). 

The AGN power is on average lowest at SFR< 0.1 M0 yr^ ^ . 
This weak AGN power in quiescent galaxies is again the point 
we have often mentioned above. Weak AGNs tend to reside in 
quiescent galaxies in our sample. In other words, AGNs on the 
red sequence are low-luminosity AGNs. These are the galaxies 
that BPT tends to miss; as we have shown in Fig. 12, the BPT 
gives a low AGN rate on the red sequence, while we find a 
large AGN population there. On the other hand, strong AGNs 
are located in actively star forming galaxies and we find that 
they tend to fill in the green valley in line with previous studies 
(Schawinski et al. 2007; Schawinski et al. 2010). 

The AGN power - SFR relation observed in Fig. 18 may 
provide a physical link between the down-sizing in galaxy star 
formation (e.g., Thomas et al. 2005) and the down-sizing in 
AGN activities (Ueda et al. 2003). More massive galaxies form 
at earlier times with shorter star formation time scales. The 
peak SFRs of galaxies are therefore higher for more massive 
galaxies. Fig 10 of Thomas et al. (2005) qualitatively illustrates 
it. If the AGN power is correlated with the SFRs as shown 
in Fig. 18 at all redshifts, then we expect to observe a peak 
of space density of stronger AGNs at higher redshifts. Ueda 
et al. (2003) actually observed that a peak of stronger AGN 
activity is located at higher redshifts. Perhaps the two down- 
sizing phenomena are linked by the co-evolution of black holes 
and host galaxies. 

To obtain deeper insights into the relationship be- 
tween the AGNs and host galaxies, we compare black 
hole growth rate (L[oh]+[oiii] agn/^^bh^ which is pro- 
portional to Mbh /Mbh), and host galaxy growth rate 

{SFR/Mstellar, wMch is M stellar /M stellar) in Fig. 19. 

L[oii]+Ioiii]agn/'^bh also serves as an indicator of the 
Eddington ratio. We observe a tight correlation be- 
tween Hq and [oil]H-[oill] (L[oii]+[oiii]AGN^'^-^LHaAGN) 
and we obtain a bolometric correction factor of 50 to 
L[oii]+[oiii]agn assuming Ltoi ~ 220LHa from Ho (2008). 
We find that most of the Oxygen-excess objects typically have 
an Eddington ratio between 10"'* and 10~^. We note that a 
bolometric correction to convert an emission line luminosity 
into bolometric luminosity is highly uncertain and probably not 
a constant correction (Ho 2008). Therefore, the Eddington ra- 
tio should not be over-interpreted. 

We are again limited by the sensitivity limit and selection 
bias, but galaxies with lower mass black holes (i.e., lower 
mass galaxies) tend to show higher black hole growth rates and 
galaxy growth rates. In any black hole mass bin, the black hole 
growth rate is low at low SFR/Mstciiar and few objects with 
high L[oii]+ioiii]/^BH are observed. The upper envelope 
of the black hole growth rate increases with increasing host 
galaxy growth rate and rapidly growing black holes are always 
hosted by rapidly growing galaxies. This clearly represents the 
co-evolution of the super-massive black holes and their host 



galaxies. Again, this is surprising because these two growths 
occur at very different physical scales. The black hole growth 
occurs on a 1 pc scale (or less), while the galaxy growth occurs 
on a 10 kpc scale. But, our result shows that black holes grow 
in growing galaxies. 

There are many recent speculations on energy feedback from 
AGNs to suppress star formation activities of the host galax- 
ies. Do we observe any hint of AGN feedback in action in 
our data? The positive correlation in the upper envelope of the 
AGN power and host galaxy SFR suggests that strong AGNs do 
not appear to suppress on-going star formation. If strong AGNs 
could sharply suppress star formation, we would not have ob- 
served the positive correlation between AGN power and host 
galaxy SFRs. However, there may be a time delay between 
AGN activities and the suppression of SFRs. Also, the appar- 
ent lack of suppressed SFRs may possibly be because we do 
not probe very strong AGNs such as QSOs, in which star for- 
mation may be suppressed (Ho 2005; Kim et al. 2006). 

On the other hand, the observed large fraction of AGNs in 
quiescent galaxies may give a support to the so-called 'radio 
mode' feedback (Croton et al. 2006; Bower et al. 2006), which 
suppresses gas cooling for further star formation and keeps red 
galaxies red. Our data, however, do not tell us whether these 
AGNs are actually suppressing the gas cooling and they have 
any direct effect on the host galaxies. We cannot make anything 
more than a speculation here, and it would be fair to say that 
we do not observe any convincing evidence of AGN feedback 
in action and we do not observe any evidence against it. 

Some of the results in this subsection are already mentioned 
in previous studies using colors as a proxy for star formation 
rate. But, we have characterized them with physical parameters 
such as host galaxy SFRs. We shall also emphasize again that 
we have studied intrinsic AGN emission thanks to the unique 
feature of the Oxygen-excess method taking into account the 
selection bias. We have put the previous results on the quan- 
titative ground and it is now clear that growing black holes 
are often hosted by growing galaxies, which represents the co- 
evolution of black holes and the host galaxies. 

5. Summary and Discussion 

In the last section, we have studied correlations between host 
galaxy properties and AGN activities. Among the numerous 
plots presented there, our primary findings can be summarized 
as follows. 

1 . The AGN fraction increases with increasing stellar mass 
of the host galaxies. Most of the AGNs are located in 
> 10^" Mq galaxies. 

2. The AGN fraction is higher in red early-type galaxies 
than in blue late-type galaxies. But, this is at least partly 
driven by the fact that the Oxygen-excess method tends 
to miss low-luminosity AGNs in actively star forming 
galaxies. 

3. Most of the identified AGNs are LINERs. The 
LINER/Seyfert partition is primarily a function of AGN 
power and its dependence on host galaxy properties is 
weak. Seyfert-like AGNs with [ Oil] ag Ar<[oni]AGAr 
appear only in powerful AGNs. 
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4. AGN power does not strongly correlate with black hole 
mass (or stellar mass of the host). 

5. Powerful AGNs are located in actively star forming 
galaxies. 

6. Rapidly growing black holes are hosted by rapidly grow- 
ing galaxies. 

7. We observe no direct evidence of the AGN feedback in 
action. 

We find the 1st and 4th points above very interesting. The 
first point is shown in Fig. 1 1 and the 4th point is in Fig. 5 and 
the bottom-right panel of Fig. 17. Mass is the primary param- 
eter to characterize the fraction of AGNs. The AGN fraction is 
very small below 10^" M©, but it increases rapidly above that 
mass. One might suspect that the lower AGN fraction at lower 
mass is because lower mass black holes show weaker activities 
and they go undetected. However, the 4th point above suggests 
that is not the case. The correlation between AGN power and 
mass is very weak and it is unlikely that such a weak correlation 
can explain the observed strong mass dependence of the AGN 
fraction. Recently, Mullaney et al. (2012) reached the same 
conclusion. Furthermore, lower-mass galaxies are located at 
lower redshifts due to the nature of a flux limited survey. The 
detection limit on L[oii]+[oiii]agn shown by the dashed lines 
in the figures is a limit at which we can detect it at all redshift 
(0.02 < z < 0.10) and our sensitivity is actually better for low- 
mass galaxies which we can observe only at low redshifts. The 
observed lack of low-mass AGNs is not due to the absence of 
bulge either As shown in Fig. 14, the AGN fraction is very 
low at low mass regardless of the morphology. While very low- 
mass AGNs do exist (e.g., NGC4395 has Mbh 4 x lO^M©; 
Peterson et al. 2005; see also the work by Barth et al. 2008), the 
observed strong mass dependence suggests that such low-mass 
AGNs are rare. It seems that mass is like a 'switch' to activate 
super-massive black holes. 

There are at least two very naive interpretations of this ab- 
sence of low-mass AGNs. One is that super-massive black 



holes exist in all galaxies, but low-mass galaxies cannot trans- 
port material to the center and feed super-massive black holes 
for some reason. That is, the majority of the black holes are 
inactive and thus are not detected. In this scenario, a process 
to trigger AGN activity is unlikely related to the presence of 
bulge because the AGN fraction is nearly zero in low-mass 
galaxies regardless of their morphology (Fig. 14). Also, the 
presence/lack of a transportation mechanism must be largely 
independent of the accretion rate as we observed no strong cor- 
relation between the AGN power and black hole mass (Figs. 5 
and 17). Namely, a mass transportation mechanism is estab- 
lished at a higher rate in more massive galaxies, but the trans- 
portation rate should be largely independent of mass. Although 
such a physical process is somewhat difficult to imagine, this 
is a possible interpretation of the result. 

The other very naive interpretation is that the majority of 
low-mass galaxies do not host super-massive black holes. This 
is perhaps an extreme idea, but it does explain the observed 
trend because the presence/absence of super-massive black 
holes works just like a switch. If a galaxy hosts a super-massive 
black hole, it can be active if the host can feed it. If a galaxy 
does not harbor a super-massive black hole, it can never be ac- 
tive regardless of the host galaxy properties. In this scenario, 
mass of the host galaxies is a probability function of the pres- 
ence of super-massive black holes. The form of the probability 
function is such that more massive galaxies are more likely to 
host super-massive black holes and the probability is very small 
at < 10^° M©. In fact, the super-massive black hole occupation 
fraction is apparently small in low-mass galaxies (e.g., M33 
and NGC205 do not host a super-massive black hole; Gebhardt 
et al. 2001; Valluri et al. 2005). This provides supportive ev- 
idence of the decreasing probability of hosting super-massive 
black holes with decreasing the host galaxy mass. 

Mullaney et al. (2012) observed that the AGN activity is 
nearly independent of host galaxy stellar mass at < z < 3. 
They also observed that AGNs are almost exclusively hosted 
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by > 10 Mq galaxies. It seems that the trends that we observe 
at z = holds out to z = 3. This might suggest that a primary 
epoch of the formation of super-massive black holes is located 
at very high redshift. Dark matter halos that collapse at high 
redshifts are among the most massive ones today and therefore 
the dominance of AGNs in massive galaxies is naturally ex- 
plained if the formation of super-massive black holes peaked 
at high redshifts. The absence of low-mass AGNs might in 
turn imply that the super-massive black hole formation is not 
very efficient at low redshift. Perhaps the mass dependence of 
the AGN fraction reflects the formation history of the super- 
massive black holes. 

Our data do not allow us to conclude which of the two sce- 
narios discussed above is more likely. Also, there may well 
be other ways to interpret our result. Before we try to make 
any further speculations, we have to make an attempt to con- 
firm that an AGN fraction in low-mass galaxies is actually very 
low compared to that in high-mass galaxies. An SDSS fiber 
subtends 3 arcsec on the sky, which corresponds to ^ 3 kpc at 
z = 0.05. AGN emission, if any, is therefore strongly contami- 
nated by host galaxies. Also, our black hole mass estimates are 
probably not very accurate. Greene & Ho (2007) argued that 
the black hole mass function decreases at low mass. Our result 
may not be directly comparable to theirs because their result is 
based on objects with broad Ha emission, while ours is based 
on narrow line AGNs. But, the observed decline in the black 
hole mass function at low-mass end is qualitatively consistent 
with ours. Still, it would be necessarily to carry out a dedicated 
observation to sample the nuclear spectra of nearby low-mass 
galaxies to achieve a better sensitivity to identify AGNs. If the 
absence of AGNs in low-mass galaxies is confirmed in such 
a survey, then we will be in a position to further discuss the 
physical origin of such a trend. 
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